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We extend and improve earlier estimates of the ability of the proposed LISA (Laser Interferometer 



Space Antenna) gravitational wave detector to place upper bounds on the graviton mass, 



by 



comparing the arrival times of gravitational and electromagnetic signals from binary star systems. 
We show that the best possible limit on m g obtainable this way is ~ 50 times better than the 
current limit set by Solar System measurements. Among currently known, well-understood binaries, 
4U1820-30 is the best for this purpose; LISA observations of 4U1820-30 should yield a limit w 3 — 4 
times better than the present Solar System bound. AM CVn-type binaries offer the prospect of 
improving the limit by a factor of 10, i/such systems can be better understood by the time of the 
LISA mission. We briefly discuss the likelihood that radio and optical searches during the next 
decade will yield binaries that more closely approach the best possible case. 
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Recent work by Will Q and Larson and Hiscock ||] 
has examined how gravitational wave observations by 
LIGO H and LISA Q can be used to place upper 
bounds on the mass of the graviton. Will's proposed 
method utilizes the dispersion of the waves generated 
in binary inspiral caused by a nonzero graviton mass, 
while Larson and Hiscock have proposed direct corre- 
lation of gravitational wave (GW) and electromagnetic 
(EM) observations of nearby white dwarf binary star 
systems. Both approaches promise improved bounds 
compared to the present bound based on Solar System 
dynamics, rrig S < 4.4 x 10~ 22 eV, which corresponds 
to a bound on the graviton Compton wavelength of 
Af 5 = h/m g = 2.8 x 10 12 km g. 

The ultimate limit on the binary star method is de- 
termined by the precision with which the phase of the 
gravitational wave signal can be measured. Larson and 
Hiscock [|| estimated this uncertainty in phase by consid- 
ering the cadence of measurements made by LISA dur- 
ing a one-year integration of the periodic signal from 
the binary star system. However, this is not the dom- 
inant source of uncertainty in the gravitational wave 
phase measurement. Let <f> be the orbital phase of the 
binary at some fiducial time. For GW measurements 
with signal-to- noise S/N >> 1, one can estimate the 
uncertainty with which this phase can be extracted by 



u GW = y/^y"t>{l + O(SfN)- 1 ), where is the 
Fisher information matrix. For a circular binary, the sig- 
nal is characterized by seven or eight physical parame- 
ters: two angles describing the direction TV to the binary 
(i.e., its position on the sky), two more angles describ- 
ing the normal L to the orbital plane, the overall am- 
plitude A, the overall phase 4>, the orbital frequency /o, 
and (non-negligible for for some binaries) the frequency 
derivative /o- If all parameters except cj) are somehow al- 
ready known (e.g. if the binary is resolved optically, for 
example by the Space Interferometry Mission (SIM) ||), 
then (r- 1 )"^ = (r^)- 1 = \{S/N)- 2 §, so we would 



7 (1 + O(SfN)- 1 ), where T l3 
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^ (jf) . More generally we can write 
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where a(> 1) is a correction term that accounts for the 
fact there will generally be additional unknown parame- 
ters to be extracted from the GW data, which can only 
increase (r -1 )*^. The case of interest, for setting a limit 
on m g , is where 6(f>aw « h or jf >> h so from now 
on we ignore the O(SfN)- 1 correction in Eq. @). 

How large is a likely to be? To answer this question, 
we must consider more carefully the definition of <f>, and 
also consider what information will be available to sup- 
plement the GW measurement. 

By </), we will mean the orbital phase, at some fiducial 
time to j measured from the point in the orbit where the 
orbital plane intersects the plane perpendicular to the 
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line-of-sight. More precisely, 



cf, = t&n-^ifito) ■ y)/(f(t ) ■ x)] 



(2) 



where f is the (unit) orbital separation vector (point- 
ing from the more massive to the less massive body), 
x = N x L/\\N x L\\, and y = L x x. Fig. |l| illustrates 
these quantities. This is a useful definition when com- 
paring to most EM measurements, where generally one 
does not know the overall orientation of the binary, but 
can measure the phase of one body's motion towards or 
away from the observer (e.g., by Doppler measurements, 
or, as discussed below with 4U1820-30, because the flux 
is greatest at the instant the the body is furthest from 
the observer). Rather than report the phase at some 
fiducial time, astronomers typically report some "epoch" 
(i.e., instant of time) to when the orbital phase (mea- 
sured optically) is 4*em = 0. Of course, this optical (or 
X-ray, etc.) measurement will have some error; we de- 
fine —S(f>EM (~ — 27r/o<5£o) to be the true EM phase at 
time to- Our proposal for constraining m g amounts to 
using LISA to measure 4>Gw{to)- If photons and gravi- 
tons travel at the same speed, then (f>Gw(to) should be 
consistent with zero. More precisely, one can set the fol- 
lowing lower bound on the Compton wavelength X g of 
the graviton (see M for details): 



(3) 



where D is the distance to the binary system and A is 
defined by p2l 




A = 



6cj> 2 



EM 
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Note that in Eq. (||) we have assumed that the errors 
5(f>EM and Stficw are uncorrelated. 

Eq. (|^) represents one convention for the zero-point 
of the phase </>; with a different convention, Eqs. (g)— (Q) 
would remain valid, but with different numerical values 
for 6<pEM and Stficw- Obviously, for any given binary, 
one wants to define <f> in a way that minimizes A. For 
instance, if the binary is seen nearly face-on (i.e., L and 
N are nearly parallel), then the <j) defined in Eq. (|4[) is 
difficult to measure, since a small shift in L can produce 
a large shift in tf>, for fixed f(to). In this case, if one 
could somehow resolve the two components of the binary 
optically (e.g., with SIM), then one would probably be 
better off choosing some arbitrary-but-easily-determined 
direction (e.g., the direction to Polaris) to define = 0. 

For the purpose of determining the best possible upper 
limit one can set on X g , we will next assume that |5<^bm| 
is small compared to |(5</>gw|- We expect this will gener- 
ally be the case when one can determine the optical phase 
at all (because of the much higher S/N one typically has 
with optical measurements) , but this will have to be ver- 
ified on a case-by-case basis. Combining Eqs. ([j]), (^), 
and (H) with this assumption then allows us to obtain a 



simple expression for the best possible limit achievable on 
the Compton wavelength of the graviton from combined 
EM and GW observations of binary systems: 



A, 
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The rms S/N (averaged over source locations and ori- 
entations) is given by fel 



S/N = h(2T)^/[S s h A {f gw )} 



1/2 



(6) 



where h = <J < h\ + h \ > is the rms value (averaged 

over time and source direction) of the strain field at the 
detector, T is the observation time, f gw — 2fo is the 
gravity-wave frequency, and Sf^ A is the "sky-averaged" 
spectral density of the detector noise. For a circular- 
orbit binary, h — (£ , gtu ) 1 / 2 /(27r/oZ?), so the numerator 
in Eq. (|J) is given by 



hV2T = 



7.52 x 10 _17 Hz" 



-1/2 
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where the stellar masses M\ and M 2 are measured in 
units of the solar mass. The sky averaged spectral density 
of detector noise is given below the transfer frequency 
/, = c/(2nL) by D 



S, 



SA 



2(S a n + S*)/R 



gw 
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where R gw = 3/5 is the low frequency gravitational wave 
transfer function for LISA, and 

8S Q 



L 2 (2nf gw y 



and 



o,. 



2S X 
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are the spectral density of acceleration and position noise 
(respectively) which are output through the detector. L 
is the armlength of the interferometer. The spectral den- 
sities S a and S x are the raw spectral noise densities of 
the noise acting on the detector. The LISA specifica- 
tions for these values are \fS^ = 3 x 10 _15 m/(s 2 v / Hz) 
and VST = 2.0 x 10- 11 m/\/Hz' §,(2|. 

Since S/N is inversely proportional to the distance to 
the binary, D, X g as given in Eq.(^|) is actually indepen- 
dent of the distance. (This independence of the limit 
on the source distance was already noted by WillQ| in a 
similar context.) It is then worthwhile to combine Eq.(|^) 
and Eq.(Q) to obtain 



1.05 x 10 4 
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Utilizing this expression for X g , and the low- frequency 
approximation to the predicted LISA sensitivity curve as 
given by Eq. (||), Q, we can now determine the best 
possible lower limit on X g that could be obtainable via 
this method, with an "ideal" binary acting as the si gna l 
source. The system variables that appear in Eq. (|ll| ) 
are the masses of the two stars, their orbital frequency, 
and the sensitivity of the gravitational wave detector 
( "noise" , to be evaluated at the frequency of the grav- 
itational wave, f gw = 2f a ). The strongest bound on X g 
will occur for a binary system whose orbital frequency 
fo is equal to the frequency f c that minimizes the func- 
tion f 2 / 3 S% A {f). Utilizing Eqs. (§0), this minimum is 
found to occur at a frequency 

f c ~ 2.06 x l(T 3 Hz , (12) 

and that the sky-averaged spectral density of detector 
noise at the corresponding f gw = 2/ c is 

S£ A {2f c ) = 1.28 x KT^Hz- 1 . (13) 

The strongest limit on X g is obtained by assuming that 
the stellar masses are equal and as large as possible; we 
will take them both to be equal to the Chandrasckhar 
mass, Mi — M% ~ 1.4M Q , which is the maximum mass 
for a white dwarf, and appears to be the "canonical" mass 
for neutron stars, based on observation. Evaluating A g 
by substituting these values into Eq. ( pi] ) , we find that for 
a 3-yr measurement, LISA could set an upper limit of 

A™ ax ~ (a~ 1/a )1.4 x 10 14 km , (14) 

which is about 50/y/a times stronger than the present 
solar system limit on X g . For more typical WD masses, 
Mi = M2 = 0.5Mq, the improvement factor is still « 
20/VS. 

Note that in Eqs. (|l2|) (p~4[) for the optimum frequency 
and optimum limit, we have included only LISA's instru- 
mental noise. This is reasonable at gravity-wave frequen- 
cies f gw > 2 mHz, but for f gw < 2 mHz, confusion noise 
due to unresolved galactic and extragalactic binaries is 
probably the dominant LISA noise source, and the total 
noise increases (sum of instrumental and confusion noise) 
rises steeply at lower frequencies. E.g., for frequencies 
2/o S 1 mHz (i.e, at GW frequencies a factor of two or 
more below where WD confusion noise begins to dom- 
inate), f^S^ (2/c)] 1 / 4 > 7/ 1/6 [S^(2/ c )] 1 /4 (where 
here we include the WD confusion noise in S^ A ). Thus 
for binaries with 2/o < 1 mHz, the best limit one could 
set on m g is only ~ better than the Solar System 

bound. For this reason, we will concentrate on binaries 
above the frequency where confusion noise dominates: 
2/o > 2 mHz. 

We now return to our discussion of a. For sources 
where an EM/GW comparison can be made, it is clear 
that the sky location N will be known to extremely high 
accuracy from the EM signal. It is reasonable to expect 
the frequency fo (at epoch to) and its derivative fo can 



also be determined from the EM data. That leaves four 
parameters, including cf>, to be determined by LISA. Fol- 
lowing the methods described in 0, we have calculated 
the Fisher matrix for this four-parameter problem. We 
find that for the ~ 20% of cases where \L ■ N\ < 0.2 (i.e., 
the cases where the binary is seen nearly edge-on), the 
degradation of 5(f>GW is quite small: a < 1.2. And for 
the best ~ 40% of cases, the degradation factor a is less 
than ~ 3. So a fair fraction of sources will be favorably 
oriented for determining 8<Pgw- Fortunately, the "edge- 
on" orientation that is favorable for small 54>qw is also 
favorable for small 6<Pem, since this orientation gives the 
largest Doppler shift. 

We next consider one source, 4U1820-30, which 
amounts to an "existence proof" of the feasibility of this 
method of constraining m g by comparing optical and 
gravitational arrival times. 



A. 4U1820-30 

Low-mass X-ray binary 4U1820-30 appears to con- 
sist of a low-mass (~ 0.07M Q ) white dwarf in orbit 
around a NS. The orbital period is 11.46 ± 0.04 min, 
so few = 2.909 mHz-i.e, a frequency where the galactic 
background can probably be subtracted out (so instru- 
mental noise dominates), and close to the optimum fre- 
quency for constraining m g . The 11.46 min periodicity 
was first observed in X-rays, but was also recently de- 
tected in the UV by the Hubble Space Telescope's Faint 
Object Spectrograph (FOS). (The angular resolution of 
HST was required for the measurement, since 4U1820-30 
is in a very crowded field, near the core of globular cluster 
NGC 6624.) The UV modulation (and roughly its ampli- 
tude) had been predicted by Arons & King Jl4] , based on 
the following picture. The WD rotation period is tidally 
locked to the orbital period, so that the same side always 
faces the NS. This is the WD's "hot side," as it is heated 
by X-rays from the NS; the UV flux we measure varies as 
the hot side is alternately facing towards and away from 
us. Clearly, the maximum UV flux occurs when we see 
the hot side most nearly straight on, which occurs at the 
point in the orbit when the NS is closest to us. This ob- 
servation provides crucial understanding of the relation 
of the phase of the binary's light curve to that of the 
associated GW signal, which is presently not understood 
for stronger and more well-known binary systems such 
as AM CVn. The measurements in Anderson et. al Jlj| 
determined the overall phase (equivalently, the epoch of 
UV maximum) to within S(/)em ~ 0.16. (They state 
<5(epoch) = 0.0002 days.) Its distance is estimated at 7.6 
kpc, which means LISA should detect it with S/N w 25 
(for a 3-yr observation, using the results from a single 
synthesized Michelson) . Based on Arons & King ]j3 , the 
binary's inclination angle i is estimated to be i w 43°±9°, 
but this is somewhat model-dependent. We have calcu- 
lated that a 1 / 2 ~ 4 for binaries with i w 45°, typically. 
Using this value for y/a in Eq. ([ll]), we estimate that 
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LISA observations of 4U1820-30 should improve the So- 
lar System bound on m g by a factor ^3 — 4. This im- 
provement is comparable to what should be obtainable 
by analysis of GW signals from inspirals of stellar-mass 
compact objects observed by ground-based interferome- 
ters such as LIGO and VIRGO (!]] . If we were to include 
constraints on the allowed range of L ■ N (from the op- 
tical measurements), when fitting to the GW data, that 
would of course decrease a and so improve the limit. 



TABLE I: Known AM CVn-type systems. 



Orbital Period yfeiX g V&^g/^g 



Name 


(s) 


(10 12 km) 




AM CVn 


1028.73 


27. 


9.7 


EC15330-1403 


1119 


23. 


8.1 


Cet3 


620.26 


25. 


9.0 


RX J1914 +24 


570 


36. 


13. 


RX J0806.3+1527 


321 


25. 


8.9 



B. AM CVn-type binaries 

Several of the "classic" AM CVn-type binary systems 
offer high potential S/N for gravitational wave observa- 
tions, as well as having sufficiently short orbital period 
to place their GW emission at frequencies high enough to 
avoid the confusion noise from Galactic and extragalac- 
tic binaries. However, these Helium cataclysmic variable 
systems, containing accretion disks, offer very compli- 
cated light curves that make it difficult to understand 
how the binary's light curve phase is related to the line 
of masses connecting the two stellar components. Un- 
less the relative phase of the EM and GW signals at the 
source is known, a binary cannot be used to place useful 
limits on the graviton mass. 

However, virtually all studies of these systems to date 
have utilized time-resolved photometry; little or no time- 
resolved spectroscopic observations have yet been dedi- 
cated to these systems. Time-resolved spectroscopic ob- 
servations should be able to provide Dopplcr information 
that will resolve the ambiguous relation between EM and 
GW phases at the source. As an example, in the epony- 
mous AM CVn system, the orbital velocity of the primary 
star is about 40 km/s; today, largely driven by the extra- 
solar planet research efforts, Doppler surveys are reaching 
accuracies of between 3—10 m/s. If such accuracies can 
be obtained in spectrographic studies of AM CVn-type 
binaries, the uncertainty in the EM phase, 8<Pem, will be 
significantly less than the uncertainty in the GW phase, 
S4>gw, for any source for which (S/N) aw i$ 1000. Since 
there is roughly a decade before LISA's launch in which 
spectrographic techniques will continue to improve, and 
such observations may be made of these binary systems, 
we feel there is a substantial chance their nature may be 
sufficiently well understood so that their EM and GW 
signals may be used to constrain X g . 

In Table 1 we display some "best limits" that might be 
obtainable from the higher frequency known AM CVn- 
type systems. The Table gives the period of the binary 
system, along with the best lower limit on X g that might 
be attainable, and the ratio of that limit to the present 
solar system bound on X g . In determining these "best 
limits" , we have assumed that optical astronomers will 
be able to adequately determine the physical elements of 
these nearby binary systems, so that the primary limita- 
tion on our method is the accuracy with which the phase 
of the gravitational wave signal can be measured. 



It is also worth noting that three of the five high- 
frequency systems listed here have only recently been 
recognized as extremely short period binaries. Cet3 (also 
known as KUV 01584-0939) was discovered in 1984, but 
its nature has been only just been revealed by hig h speed 
photometry 0. Similarly for RX J1914+24 @ £§ and 
RX J0806.3+1527 @. This suggests that additional 
such systems may well be discovered (or recognized as 
such) before LISA's launch. 



C. Prospects for discovering a binary pulsar with 

P b < 1000 s 

The discovery of a pulsar in a short-period (P b < 1000 
s) binary with a NS or WD companion would likely pro- 
vide an excellent system for constraining m g , for two rea- 
sons. First, a higher-mass system tends to give a stronger 
limit on m g . Second, relativistic corrections to the binary 
orbit (perihelion precession and orbital inspiral) and the 
pulse arrival time (the Einstein delay and Shapiro de- 
lay), often allow most of the binary's parameters to be 
extracted. In the best cases, all binary parameters are 
extracted, except for one angle: the direction of L x N. 
So only two parameters, 4>cw and this direction angle, 
need to be extracted from the GW data, which should 
generally translate into a low a. 

The two shortest-period binary radio-pulsar systems 
currently known (and also having companion mass M 2 > 
0.1M©) are are J0024-72W (P b = 2.6hr, M 2 = 0.15M Q , 
d = 4.5 kpc) and B1744-24A (P b = 1.9hr, M 2 = O.LM , 
d = 7.1 kpc)- within a factor 16 and 12, respectively, 
of our ideal period l// c = 0.16 h. (See Table 4 in 
Lorimer Q.) Recently discovered PSR J1141-6545 is 
also notable in this context, because in addition to having 
a short-period (P b — 4.74 hr), the mass of the compan- 
ion WD is rather large: M 2 > 1.0M Q J2(|. Yungelson et 
al. ]l7j] estimate that our Galactic disk contains several 
tens of NS's in binary systems with fo > 1 mHz (~ 40 
NS-WD's and - 10 NS-NS's), so such short-period bi- 
nary pulsars are likely to exist. Until now there has been 
a strong selection effect against finding short-period bi- 
nary pulsars, since the orbital motion smears out the 
pulse frequency, and the "acceleration searches" tradi- 
tionally used to demodulate the signal are ill-suited to 
observations lasting longer than ~ P b /2n. Significantly 
more sophisticated search strategies are now being implc- 
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mented (see Jouteux et al. Q and references therein), 
so it is reasonable to expect a significant extension to our 
database of binary pulsars, at the short-period end. 

Even more promising are NS binaries in globular clus- 
ters. Benacquista, Zwart, & Rasio [|l] estimate that 
Galactic globular clusters contain ~ 2 NS-NS and ~ 10 
NS-WD binaries (with a factor ~ 10 uncertainty in ei- 
ther direction) with /o > 1 mHz. Once LISA has de- 
tected these systems, LISA's few-degree angular resolu- 
tion should allow the host globular cluster to be identified 
uniquely And since LISA will provide the orbital period 
and phase to high precision, there will be only one or- 
bital parameter to search over (the maximum velocity of 
the NS along the line of sight), greatly facilitating radio 
identification of any such sources that are first discovered 
by LISA. 

D. Conclusions 

We have shown that correlating EM and GW (LISA) 
observations of LMXB 4U1820-30 should improve the 



current Solar System bound on m g by a factor ^3 — 4. 
We showed that for an "ideal" source, the improvement 
factor would be ~ 50. Since the bound on m g that one 
obtains is independent of the distance to the source, it 
seems almost inevitable that future EM surveys with in- 
creased sensitivity will reveal new (generally more dis- 
tant) binaries that more closely approach the ideal im- 
provement factor. 
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FIG. 1: The coordinate system used to define the orbital 
phase is determined by the line of sight to the binary N and 
the binary orbital angular momentum L, such that x points 
along the line of nodes, and y lies in the orbital plane. The po- 
sition vector f points along the binary axis, from the primary 
(largest mass , Mi) to the secondary (smallest mass, M2) in 
the system. The orbital phase cj> is defined to be the angle in 
the orbital plane between the +x axis and the position vector. 



